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Abstract 

We classify (l+3)-dimensional Schrodinger equations for a parti- 
cle interacting with the electromagnetic field that are solvable by the 
method of separation of variables. As a result, we get eleven classes 
of the electromagnetic vector potentials of the electromagnetic field 
A{t,x) = {AQ{t,x), A(t,x)) providing separability of the correspond- 
ing Schrodinger equations. It is established, in particular, that the 
necessary condition for the Schrodinger equation to be separable is 
that the magnetic field must be independent of the spatial variables. 
Next, we prove that any Schrodinger equation admitting variable sepa- 
ration into second-order ordinary differential equations can be reduced 
to one of the eleven separable Schrodinger equations mentioned above 
and carry out variable separation in the latter. Furthermore, we apply 
the results obtained for separating variables in the Hamilton-Jacobi 
equation. 

I. Introduction 

Being invented by Fourier and Euler long ago, the method of separation of 
variables is still the most powerful and efficient one for integrating linear 
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partial differential equations (PDEs). This is especially the case for PDEs 
having variable coefficients, where the standard Fourier transform is no longer 
applicable. Moreover, this method proves to be a useful tool for constructing 
particular solutions of some nonlinear partial differential equations such as 
the nonlinear Laplace [^, wave 0, H, and heat conductivity equations 
in (1+1) dimensions. 

The principal object of study in the present paper is a problem of the sep- 
aration of variables in the Schrodinger equation (SE) for a particle interacting 
with the electromagnetic field, 

ipO-PaPa)i'it,x) =0. (1) 

Here we use the notations 

d d 
Po = iw: - eAo(t,x), Pa = i^ eAa{t,x), a =1,2,3, 

Ot OXa 

where A = {Aq, Ai, A2, A^) is the vector potential of the electromagnetic 
field, e = const. Hereafter the summation over the repeated indices from 1 
to 3 is understood. 

Bocher is believed to be the first to obtain in [0 a systematic classifi- 
cation of coordinate systems enabling separability of the three-dimensional 
stationary SE, 

(-A + E)V^(f) = 0, (2) 

where A is the Laplacian in three dimensions. 

He has shown that equation is separable via the separation Ansatz, 

ipix) = (pi{uJi{x))(p2iuJ2ix))(p3{uJ3{x)) (3) 

in eleven inequivalent coordinate systems uJi{x), u!2{x), uj^^x). We list these 
following the famous paper by Eisenhart [§, where a rigorous geometric 
derivation of the corresponding results is given. Note that the coordinate 
systems are given in the implicit form Xa = Za{u!i, uj2,uj3), a = 1, 2, 3. 

1. Cartesian coordinate system, 

Zl = UJi, Z2 = i02, Z3 = UJ3, 
UJl,UJ2,UJ3 G R. 

2. Cylindrical coordinate system. 
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zi = e"^^ C0SCJ2, Z2 = e'^^ sina;2, Z3 = u^, 
< U2 < 2-71, Ui, UJ3 G R. 

3. Parabolic cylindrical coordinate system, 

Zi = (C^i - UjI)/2, Z2 = UJlLLt2, Z3 = UJ3, 

> 0, UJ2, G R. 

4. Elliptic cylindrical coordinate system, 

;2i = a coshtc;i costc;2, ,22 = asinhwi sinci;2, z^ = uj^, 
LO\ > 0, — TT < a;2 < TT, Go's G R, a > 0. 

5. Spherical coordinate system, 
Zi — uji^ sech U2 cos a;3, 

-^2 = <^r^ sech a;2 sin us, 
Zs — cui^ tanha;2, 

a;i > 0, a;2 G R, < a;3 < 27r. 

6. Prolate spheroidal coordinate system, 
Zi — a csch cui sech 002 cos cu^, a > 0, 
Z2 — a csch a;i sech 002 sin a;3, 

2^3 = acotha;itanha;2, (4) 
a;i > 0, a;2 G R, < a;3 < 27r. 

7. Oblate spheroidal coordinate system, 
zi — a CSC a;i sech CU2 cos a;3, a > 0, 
Z2 — a CSC a;i sech a;2 sin uj^, 

Zs — a cot a;i tanh 002, 

< a;i < 7r/2, a;2 G R, < a;3 < 27r. 

8. Parabolic coordinate system, 

= e'^^ ^'^^ cos Us, Z2^ e'^i+'^2 ^ 
^3 = (e2-i - e2-^)/2, 

ciJi, u;2 £ R, < < 27r. 

9. Paraboloidal coordinate system, 

zi = 2acosha;i cosa;2 sinhcus, a > 0, 
Z2 — 2asinha;isina;2Cosha;3, 
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= a(cosh2c(Ji + cos2ct;2 — cosh2c(j3)/2, 
cuijCUs G R, < ct;2 < vr. 

10. Ellipsoidal coordinate system, 

Zi = a — 7 — - dn(u)2, k') sn( uj^, k), a > 0, 

sn{uji, k) 

dn(u;i, A;) 

Z2 = a — — cn(u;2, k ) cn(u;3, k), 

sn(c<Ji, k) 

cn(a;i, k) 

Z3 = a — — sn(u;2, k ) dn(u;3, k), 

sn^tui, k) 

0<uJi<K, -K' <UJ2<K', < cjs < AK. 

11. Conical coordinate system, 
Zi = uJi^dn(uj2, k') sn(co'3, /c), 
Z2 = uj^^cn{uj2, A;')cn(co'3, k), 
Z3 = uj^^sn{uj2, k') dn(co'3, fc), 

ui > 0, -K' <UJ2<K', < CJ3 < AK. 

Here we use the usual notations for the trigonometric, hyperbolic and Jacobi 
elliptic functions, k{0 < k < 1) being the modulus of the latter and k' = 

(1 - fc2)V2. 

By the evident reasons, the coordinate systems 1, 2-4, and 5-11 from the 
above list are called completely split, partially split, and non-split, corre- 
spondingly. 

Note that the above list differs slightly from the one presented in since 
we have rearranged the coordinate systems u!i,uj2,uj3 in such a way that the 
relations 

Aa;„ = 0, a =1,2, 3 (5) 

hold for all the cases 1-11. 

Next, for each of the coordinate systems, Eisenhart determined the 
form of the potential V{x) that permits the separation of variables. 

The Eisenhart's technique has been applied by Olevskii [|T^] for sepa- 
rating variables in the Laplace-Beltrami operator in the spaces of constant 



curvature (see, also 11 



Smorodinsky and Winternitz with co-workers |T2|, |13[ started a systematic 
study of potentials for which the stationary SE in two and three dimensions 
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admits the separation of variables in two or more coordinate systems (so- 
called superintegrable potentials). The classification of these potentials has 
been completed by Evans 

In the mid-seventies a series of papers by Miller and Kalnins appears, 
where a symmetry approach to variable separation has been developed. This 
approach is based on the well-known fact that a solution with separated 
variables is a common eigenfunction of first- or second-order differential op- 
erators, which commute with each other and with the operator of an equa- 
tion under consideration. Further details and an extensive list of references 
can be found in the monograph (also see, the review by Koornwinder 
Boyer, Kalnins, and Miller have obtained a systematic solution of the 
variable separation problem in the time-dependent (H-2)-dimensional free 
SE within the framework of the symmetry approach Later on, Boyer 



has described all coordinate systems providing separability of the (1+2)- 
dimensional SE having the potential V{xi,X2) = a/x\ + [3/x\ |T^. Reid has 
completely solved the problem of variable separation in the three-dimensional 
time-dependent SE for a free particle [0. 

Independently, the symmetry approach to the separation of variables in 
the equations of quantum mechanics and quantum field theory was devel- 
oped by Shapovalov, who was the first to give a systematic treatment of 
the problem of variable separation in the Dirac equation using its non-Lie 
symmetry and by Bagrov with collaborators (see, and references 
therein). Shapovalov and Sukhomlin have obtained some separable SEs 
of the form (|l|), however, their results are not complete. Let us also mention 
the papers [|^, where the physical aspects of the problem of separation 
of variables in some (1+3)- and (l+l)-dimensional SE with time dependent 
potentials are studied. 

In the present paper we are mainly devoted to the inverse problem of 
variable separation in SE (P, namely, to one of the classifying PDEs of 
the form (|I]) that can be solved by the method of separation of variables. 
Clearly, to be able to handle this problem efficiently, we need a precise al- 
gorithmic definition of what the separation of variables is. We have sug- 
gested a possible definition of the separation of variables applicable both 
to linear |2^, 2^ and nonlinear [27| PDEs, which enables developing an 



efficient approach to solving classification problems for (l+2)-dimensional 
time-dependent Schrodinger equations having time-independent [^6] , [27[| and 
time-dependent scalar potentials [p3]. Recently, we have obtained an exhaus- 
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tive classification of separable Schrodinger equations (|1]) in (1+2) dimensions 
with A3 = |2l. 



In the present paper we solve completely the problem of variable sepa- 
ration in SE dH) into second-order ordinary differential equations in a sense 
that we obtain all possible forms of the vector potential A{t,x) = {AQ{t,x), 
Ai{t,x), A2{t,x), A^ltjx)) providing separability of ([l|). Furthermore, we 
construct inequivalent coordinate systems enabling us to separate variables 
in the corresponding SEs and carry out variable separation. 

The paper has the following structure. In the second section we solve the 
classification problem for SE ([1|) and obtain all vector functions A{t,x) such 
that (|I|) is separable. What it is more, we consider briefly the problem of 
the separation of variables in SE having the generalized Coulomb potential 
and construct all the possible coordinate systems providing its separability. 
The third Section is devoted to the application of the results of Section II 
to separate variables in the (l+3)-dimensional Hamilton- Jacobi equation. 
In a concluding section we indicate some further applications of the results 
obtained in the paper. 



II. Separation of variables in the Schrodinger 
equation 

With all the variety of approaches to a separation of variables in PDEs, one 
can notice the three generic principles respected by all of them, namely, the 
following 

1. Representation of a solution to be found in a separated (factorized) 
form via several functions of one variable. 

2. Requirement that the above mentioned functions of one variable should 
satisfy some ordinary differential equations. 

3. Dependence of the so found solution on several arbitrary (continuous 
or discrete) parameters, called spectral parameters or separation con- 
stants. 



By a proper formalizing the above features we have formulated in 
an algorithm for variable separation in linear PDEs. Below we apply this 
algorithm in order to classify separable SEs of the form (|I]). 
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We say that SE (|1]) is separable in a coordinate system t, uJa = u!a(t,x), 
a = 1, 2, 3 if the separation Ansatz, 

3 

Ipit, X) = Q{t, x)(fo{t) Yl fa {uJa{t, x), X) , (6) 

a=l 

reduces PDE (|I]) to four ordinary differential equations for the functions 

iff,, (/i = 0,1,2,3), 

fo = Uo{t, ifo; A), if" = Ua{uJa, (fia, 'f'a\ (7) 

Here [/q, . . . , f/3 are some smooth functions of the indicated variables, A = 
(Ai, A2, A3) G A = {an open domain in R^} are separation constants and, 
what is more. 



rank 





dUo 


dUo 


d\i 


d\2 


dXs 


dUi 


dUi 


dUi 


d\i 


d\2 




dU2 


dU2 


dU2 


dXi 


d\2 




dUs 


dUz 


dU3 


dXi 


d\2 





The above condition secures the essential dependence of a solution with sep- 
arated variables on the separation constants A. 

Formulas (^-(§) form the input data of the method. The principal steps 
of the procedure of variable separation in SE (0) are as follows. 

1. We insert the Ansatz (|) into SE and express the derivatives ipQ, ip", 
iP2, in terms of functions (po, (pi, (p2, v^s, (p'l, ^'2-, V^s using equations 

2. We regard yjo; V^i? ^2-, V'a, V^'i, f'21 V^s^ -^i' -^25 A3 as the new indepen- 
dent variables yi,. . ., yio- As the functions Q, Ui, U2, cus are indepen- 
dent of the variables . . . , ?/io we can split by these variables and get 
an overdetermined system of nonlinear partial differential equations for 
unknown functions Q, Ui, UJ2, u;^. 

3. After solving the above system we get an exhaustive description of 
coordinate systems providing separability of SE. 
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Having performed the first two steps of the above algorithm, we arrive at 
the conclusion that the separation equations (|^) are linear both in <fo, . . . , <^3 
and Ai, A2, A3 (the principal reason for this is the fact that SE (|I]) is linear). 

Next, we introduce an equivalence relation S on the set of all coordinate 
systems providing the separability of SE. We say that two coordinate systems 
t,uji,uj2,uj3 and i,uji,LJ2,LJ3 are equivalent if the corresponding Ansatzes (|^) 
are transformed one into another by 

• the group transformations from the Lie transformation group admitted 
by SE (P with A = 0, 

• the invertible transformations of the form 

t^i = fo{t), Ui^oJi = fi{uji), (9) 
Q^Q = Qh{t)liMl2MhM, (10) 

where fo, ■ ■ ■ , fs, Iq, . . . ,1^ are some smooth functions and 2 = 1,2, 3. 

This equivalence relation reflects the freedom in the choice of the func- 
tions Q,uJi,uj2,uj3 and separation constants Ai, A2, A3 preserving the form of 
the separation Ansatz (^. It splits the set of all possible coordinate systems 
into equivalence classes. In a sequel, when presenting the lists of coordi- 
nate systems enabling us to separate variables in SE we will give only one 
representative for each equivalence class. 

Within the equivalence relation S we can always choose the reduced equa- 
tions (0) to be 

iip'o = (To(t) - Ti(t)\i) ipo, (f'^ = (Fao(cUa) + Fai{uJa)\i) (fa, (H) 

where TQ,Ti, Fao, Fai are some smooth functions of the indicated variables, 
a = 1,2,3. With this remark the system of nonlinear PDEs for unknown 
functions Q, Ui, U2, uj^ takes the form 

(12) 
(13) 

0, a = 1,2, 3; (14) 



-^-^ = 0, 1^3, ^,J = 1,2,3; 

OXa OXa 



i=i ^"^3 ^"^i 



^dxj 



+ ieQAj 



Tait), a = 1,2, 3; 

duJa ^ ( ■ du!a a ^ 

oxi \ at 
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QZ^ FioiuJi) ^ ^ + t— + AQ + 2teAa — 

~^ OXj OXj at OXa 

+g (^To(t) + - eAo - e^A.A,^ = 0. (15) 

Thus, the problem of variable separation in SE reduces to integrating the 
system of ten nonlinear PDEs for four functions. What is more, some coef- 
ficients are arbitrary functions that should be determined while integrating 
equations (p!^)- ([T5|) . We have succeeded in constructing the general solu- 
tion of the later, which yields, in particular, all possible vector potentials 
A{t, x) = {Ao(t, x), . . . , As{t,x)) such that SE (|I]) is solvable by the method 
of the separation of variables. 

The integration procedure relies heavily upon the results on the separation 
of variables in the stationary SE @. That is why we will briefly consider 
the principal steps of application of our approach for separating variables in 
equation @ with E = 1, 

AgV' - ^ = 0. (16) 
Inserting the separation Ansatz, 

lP(x) = Q{x)Lpi{uJi{x))Lp2iuJ2ix))Lp3{uJ3{x)), 

into (|T^) and taking into account the relations 

V^a = (^al(t^a) + i^a2(t^a)Al + Fa3(^^a)A2) ^a, a =1,2, 3, 

we get the system of nonlinear partial differential equations for the functions 

Q, Uli, UJ2, UJ3, 

i^) 1^1^ = 0, ^^J, .,J = 1,2,3; 

OXa OXa 

duji duJi 

[ii) 2^ Fia[uJi) — — = 0, a = 2, 3; 
t^i OXj dxj 

[lii) 2^^ + QAuJa = 0, a = 1,2,3; (17) 

i=l OXj OXj 
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Now we can utilize the classical results on variable separation in the 
stationary SE ([TBD . According to the general solution uj = uj{x) of system 
( PTD splits into eleven inequivalent classes whose representatives are given in 



As Ui,uj2,uj3 are functionally independent, the inequality 

3 



det 



duji 



dxa 



7^0 



(18) 



i,a=l 



holds. Taking into account this fact and relations (^), we get from Hi that 
Q{x) = const, whence it follows that without losing generality we may choose 
Q{x) = 1. In view of this we rewrite system (^) in the following way: 



duJi du 



' = 



= 6 



1,2,3; 



1,2,3, 



(19) 



i=l 



where 6ij is the Kronecker symbol. 

System ([T9| ) coincides with the system of equations (|12D and ([13|) under 
Ti = 1, T2 = T3 = 0. Next, in view of arbitrariness of the choice of the sep- 
aration constants Ai,A2,A3 we can always suppose that Ti(to) = l,7'i(io) = 
^2(^0) = for some to e R. Consequently, system ([T^) , (pT^) with t = to 
takes the form (|TD|). This is a key point enabling us to use the results of 
for integrating system (p!2D, (p!3D. 



Lemma 1 The general solution uo = u{t,x) of the system of partial differ- 
ential equations ^^), ( [7^ j is given implicitly by the following formulas: 



x = T{t)H{t)z{Q)+w{t). 
Here T{t) is the time- dependent 3x3 orthogonal matrix: 



I cos a cos /5 — sin a sin [3 cos 7 
sin a cos (3 + cos a sin (3 cos 7 
\ sin (3 sin 7 

— cos a sin /5 — sin a cos (3 cos 7 
-> — sin a sin /5 + cos a cos /5 cos 7 

cos/5 sin 7 



(20) 



(21) 



sma sm7 
- cos a sin 7 
cos 7 
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C(,f3,'y being arbitrary smooth functions oft; z = z{u) is given by one of the 
eleven formulas from H{t) is the 3x3 diagonal matrix, 

I hit) \ 

Hit) = h^it) , (22) 
V h,it) J 

where 

(a) hiit), h2it), h^it) are arbitrary smooth functions for the completely split 
coordinate system (case 1 from ^)), 

(b) hiit) = h2it), hiit), h^it) being arbitrary smooth functions, for the par- 
tially split coordinate systems (cases 2-4 from 

(c) hiit) = h2it) = h^it), hiit) being an arbitrary smooth function, for 
non-split coordinate systems (cases 5-11 from (^) 

and wit) stands for the vector-column whose entries wiit),W2it),W3it) are 
arbitrary smooth functions oft. 

Proof. First we perform the hodograph transformation in the system of 

PDEs riT2i), m 



t = t, Xa = Uait,UJi,UJ2,UJ3), = 1,2,3. (23) 

Direct computation shows that the following identities hold: 

1. ujix^ujj^^ = ^i^ik^ljk - ^ij^kk), ii,j, k) = cycle (1, 2, 3), 

2. uji^^uji^^ = ^iQjjQkk - ^%), ii,j, k) = cycle (1, 2, 3), (24) 



52 
where 



ilij = ^— ^— , c) = det 



dUa 



dujb 



3 

7^0. 



a,b=l 

Consequently, the initial system (p!2D , (|13|) , after being rewritten in the 
new variables, reads as 

^^ = 0, t^j, ^,^ = 1,2,3; 
E F,M)^^^^ = S^TJt), a =1,2,3. 

f ■ ,\ 1 /I o Q^ (y^i UUO^ OUk OUOk 

(j,j,fc)=cycle (1,2,3) J 7 
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If we introduce three vectors, 



du 
duji 



V2 



du 



V3 



du 



then the first three equations of system ( psD read as VaVb = 0, (a, b = 
1,2, 3, a 7^ b). Consequently, vi,V2,V3 form an orthogonal system of vec- 
tors in the space R^. It is well known from the analytical geometry that the 
most general form of these vectors can be expressed via the Euler angles, 



du 

duJi 

du 
duj2 

du 
duj^ 



cos /i cos /s - sin /i sin /s cos /a 
Ri { sin /i cos /2 + cos /i sin cos /s 
sin /a sin /a 

- cos /i sin /2 - sin /i cos /2 cos /a 
R^l - sin /i sin /2 + cos /i cos /2 cos /s 
cos /2 sin /s 

sin /i sin /g 
{ - cos /i sin /s 
cos /a 



(26) 



where /i, . . . ,i?3 are arbitrary smooth functions of t, cui, ci;2, cua- The above 
formulas give the most general form of functions dua/duh, a, 6 = 1, 2, 3 that 
satisfy the first three equations from (|2F 
remaining equations from (^) yields 



Next, inserting (|26|) into the 



Y.Fij{ui)Rf = T,{t), J = 1,2,3. 



(27) 



Thus, we have transformed system 



to equivalent form 



where one should take into account the compatibility conditions for the sys- 
tem of PDEs (HD, 

d ( dui \ 



d ( dm 



dujk \duJj / dujj \dujk 



j^k, i,j,k 



1,2,3. 



Hence, we get the overdetermined system of nonlinear PDEs for the functions 

fl, /2, /s; 



, dfi ^ dh 
duJi diOi 



-R. 



^idR, 
dujo 
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COS tr> 


+ 9f2 _ 
duj2 


^ ^dR2 
= R7 




COS /q 


+ 9f2 _ 
dujs 


= 0, 






COS To 

OUJi 


+ sin f-y 


sm Tr> — 





r dh 

COS To 


+ sin fo 


. , dh 

sm = 

du;2 


— -TLq 

' 9u;3 


COS h r. 


+ sin /2 


. , dh 

sm/3 = 
da;3 


^-l9i?3 


. f dh 

ouJi 


- cos h 


. rdh 

sm /a = 


^.^dR, 
' dus 


. r dh 

sm h 
c'a;2 


- COS h 


. fdh 

sm /a = 

0UJ2 


0, 


sm h a 

OUJs 


- COS /2 


. .dh 
sm /a = 

dujs 


^-l9i?3 



While integrating the above system we have to differentiate between the 
two cases: (1) sin /a 7^ and (2) sin/3 = 0. 

Case 1. Suppose the condition sin /a 7^ holds. Then we can solve the 
above system with respect to dfa/dujh, (a, b = 1,2, 3) and get 

^ = -R3 ^— cos h CSC /a, 
dh j,-idR2 . „ „ 

-K— = --^3 ^ Sm/2CSc/3, 

aa;2 a6i;a 

dfi I dRs . 1 9i?a „ „ 

-— = i?2 ^ — sm /2 CSC /a - - — cos /2 esc /a, 
aa;a ac(;2 c/o;! 

dh j,-idR, ,dR, 

= -R2 ^ ^3 ^— cos /2 cot /a, 

dooi duj2 doos 

^ = -Ri ^ h-Ra ^— sm/2Cot/a, (28) 

-— = -i?2 ^ sm /2 cot /3 + cos /2 cot /3, 

aa;3 d(jj2 dui 

dh j,-idR, . 

13 



df3 dR 



-^s'i^ COS/2, 



duj2 duj: 



dfs jy-i dRs , , „_i ORs . 
^— = R2 ^ — COS/2 + i?i ^ — sm/2. 

OUJ3 OUJ2 OUOi 

From the compatibility conditions of the above system of PDEs we get 
the system of nonhnear differential equations for Ri, R2, R3: 

fj2 p 

1. nin2-x — 7; -tt-i — 7; -K2 — 7; — — u, 

OUJ1OUJ2 

2. R,R,^^-R,1ZZ±112L^R^:^_Z±=0 

OUJ20UJ3 

<J- -Ki-KsT^ — -rts — -rti — — — — U 

acjiacus 0UJ3 ouoi " ' " ' 

4 ^2 ^ ^ ^ ^2 ^2 ^ ^ _ ^2^ ^2 ^ ^ ^29) 

^ ^ Coco's ScJs ^ ^ 9(^2 <9co'2 ^ ^ duJi duji 

p2 p2 p 9^R2 p2 p p2 d^R-i _ n 

— it]^ /t2 -Ks 7- — K2 rc^ -X — 7^ — — u, 

OUsOUJs OUJ2OUJ2 

5 ^2 _ ^2 j^2 _^ j:^2 j:^2 ^-^3 

p2 p2 p p p2 p2 _ n 

— It^ it2 -Ks T. ^ -Kl -K2 -r<3 ^ H " U, 



dR2 


dRs 


duJi 


duj2 


dRi 


dR3 


dujs 


duj2 


dRi 


dR2 



dRi 


dR3 


duj2 


duJi 


dRi 


dR2 


duj2 


dujs 


dR2 


dRs 


duj3 





^2 j^2 _^ j^2 j^2 _^ j^2 j^2 9R2 



— i?^ i?2 -R3 7^ — Ri RI -R3 — TT— — 0. 

OUJ2OUJ2 ooJiOuJi 

Remarkably, there no need for a direct integrating of the system of nonlin- 
ear PDEs (p9D, since it has been solved by Eisenhart [§] under the assumption 
that Ri,R2,R3 are independent of t. The only thing to be done is to find 
out in which way the temporal variable t enters into the general solution of 
system (^)-(|29|) in the case under study. 

As shown above, the problem of variable separation in SE reduces to 
integrating system (|19D and, what is more, the general solution of the latter 
is given within the equivalence relation (P) by one of the formulas from @). 
Consequently, if we fix the temporal variable t to be equal to to £ then 
after integrating equations (p^)-(|29|) we get within the equivalence relation 
(|9D formulas M). 
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In view of this fact we can solve relations (|25| ) under t = to with respect 
to Fij{uji) (note that Fij are independent of t) for each class of functions 
X = z{uj) given in (^). The results of these calculations are presented below 
in the form of 3 x 3 Stackel matrices J-'i, . . . , J^n, whose (i, j)th entry is the 



corresponding function Fj, 



[UJi 



We give the canonical forms of the matrices 



J-'i, ... , T\\ up to the choice of separation constants Aj, i = 1,2,3 in (pi]). 



^6 




1 




/ G? COsh^ i)J\ 1 



-a^cos^ci;2 — 1 




Co"! 








cosh" 




^2 



sinh ^ tui 

UJ2 



a? cosh ^ 







— sinh uoi 

COSh^^ UJ2 





a? sin ^ tui — sin ^ cji 



— cosh ijj2 




cosh a;2 







3 4a;i 





2wi 



-e 




cosh^ 2u;i —a cosh 2ci;i 
— cos^ 2uj2 a cos 2u;2 
cosh^ 2(jj7i a cosh 2a;.ci 




(30) 
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10 



sn'^ , fc) sn^ , A;) 

V cn^(co'3, fc) A;^ cn^(c<;3, A;) 



1 / 



^1 



11 



^1 






-^1 
k'^cn^{uj2,k') 

A;^cn^(co'3, A;) 







With the exphcit forms of the functions F^, {i, j = 1, 2, 3) in hand, we can 
solve (^) with respect to Ri, R2, R3. Inserting the result obtained into (^) 
and splitting by c<Ji, 072, 003 yield the final forms of the functions Ri, R2, R3: 



1. 
2. 
3. 
4. 
5. 
6. 



Rl 

Rl 

Rl 
Rl 

R2 
Rl 

7. Rl 

r\ 
Rl 

8. Rl 

Rl 

^\ 
Rl 



9. 



i?9 



Rl 
r\ 



= 1,2,3; 

Rl- 

-1„2 



y-i. 



-1. 



a (cosh 2071 — cos 2c(j2) 



R3 - ^3 

T^^uj^'^, R^ = R^ = uj^"^ cosh~'^ UJ2] 
T^^a^ sinh^^ Ci;i(sinh~^ cui + cosh~^ iU2), 
T{^o? cosh"^ c<j2(sinh~^ uji + cosh~^ 102) 
T^^a^ sinh~^ uji cosh~^ UJ2] 

12 2 / 2 —2 \ 

a sin^ ci;i(sin^ cji — cosh uj2)., 



cuil^sm 

'1„2 „„„i,-2 



a cosh ^C(j2(sin ^cji— cosh 
T^^a^ sin~^ cui cosh~^ c<j2; 

Tj^^'^a^ (cosh 2c<Ji — cos 2co'2) (cosh 2c(Ji + cosh2co'3 
Tj^^'^a^ (cosh 2ci;i — cos 20^2) (cos 2u;2 + cosh 2^73), 
T^^a^ (cosh 2lji + cosh 2ti;3) (cos 2ci;2 + cosh 2^173 



10. Ri = Ti 



-1 



'dn^(co'i, k) 
sn2(co'i, k) 
dn^(co'i, k) 
sn^(u;i, k) 



-k'^cn\uJ2,k') 



'dn^(co'i. A;) 
^sn2(t^i, k) 



(31) 



+ A;^cn^(ti;3, k) 



- k"cn\u2, k') {k"cn\u2, k') + k'cn\uj3, k) 
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Rl 



Tr' {^^^^ + kVioo,, k)^ {k'\n\uj„ k') + k\n\uj,, k) 



sn^(ti;i, k) 

11. Rl = Tf ^cjf ^ Rl = Rl = T{^ujf (k''^cn^{u2, k') + k'^CYi^{u^, k)] 

Inserting the above expressions into (^), we see that the system obtained 
(we denote it temporarily as S) does not depend exphcitly on t. It is in 
involution and hence its general solution depends on three arbitrary functions 
of t. As a direct check shows the functions 

/i = arccot (^cos7 cot(/i + /?) + sin7 cot /a csc(/i + + a, 

/a = arccot (cos fs cot(/i + (3) + sin /g cot 7 csc(/i + /?)) + /a, (32) 

/s = arccos (cos 7 cos /s - sin 7 sin /a cos(/i + /5)) , 

where a, P, 7 are arbitrary functions of t and /i, /a, /s are (time-independent) 
solutions of the system S under t = t^, satisfy the system S identically. Since 
formulas (^) contain three arbitrary functions of t, they give the general 
solution of the system S. Substituting (R2|) into (B^) yields 



du 
duji 



du 
diOo 



du 

du-. 



RiT{t) 



R^Tit) 



R^T{t) 



cos /i COS /2 - sin /i sin /s cos /a 
sin /i cos /2 + cos /i ^sin /a cos /g 
sin /a sin /g 

- cos /i sin /2 - sin /i cos /a cos /s 

- sin /i sin /a +^cos /i^ cos /a cos /s 

cos /a sin /a 

sin /i sin /g 



(33) 



- cos /i sin /s 
cos/3 

where the matrix T(t) is given by formula (|21|). 

If we choose in the system obtained t = to, then its general solution is 
given within the equivalence relation S (P) by one of the formulas from (^). 
In view of this fact it is not difficult to integrate system ( ^3]) and thus get 



formulas (|20|), where hi{t), h2{t), hs{t) are expressed via Ti{t),T2{t),Ts{t), 

1. Ti = h-\ 



-4. 
11. 



T2 
T2 



1,2,3; 
= 0, 

= T3 



T3- 
= 0. 



(34) 



17 



Case 2. Let the relation sin /s = be valid. In this case we have an analog 
of system (pSD , 

dg _ ^-idRi dg__ ,dR2 



duji duj2 duj2 dcui 

dg dR, dR2 „ dR, dRs 

aciJa ouj-i ouo^ oooi 0002 

with g = ±/i + /2 and an analog of system (|29|), 

,2 d2 ^-^1 ^-^2 , n2 d2 ^-^1 "9-^2 „2 t-. t-.2 ^^^-^1 



ri]^ 1X3 -^r ;r h 1X2 iX3 — :r — JXi 1X2 -K; 



duj2 duj2 did I did I duj2duj2 

-^i^2^3^^ = 0' (36) 



atus 0UJ3 ouji 0UJ2 



System of PDEs (|35|) , (|36|) is fairly simple and is easily integrated. As a 
result we get a particular case of ( pO]) with sin 7 = 0. The lemma is proved. 

> 

With this result in hand it is not difficult to integrate the remaining 
equations from the system under study. Indeed, equations ([T^) and (p!5|) 
may be treated as algebraic equations for the functions Aj{t,x), j = 1,2,3 
and AQ{t,x), correspondingly. 

There are two different configurations of the electromagnetic field that 
should be considered separately. The first one is the case of a vanishing 
magnetic field H = rot A, namely, the case when 

Ai:,. = Aj^^, i ^ j, i,j = 1, 2, 3. (37) 

Provided the above equality does not hold we have a nonvanishing mag- 
netic field H. 

1. The case of nonvanishing magnetic field. 

Let us represent the complex- valued function Q in (||) as Q = exp(5'i -|- 
182), where 81,82 are real- valued functions. Now, if we take into account 
that the components of the vector potential A{t,x) and functions ui,uj2,uj3 
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are real- valued functions, then after inserting Q into (|14|) with the use of (|^) 
we can split the obtained equations into real and imaginary parts: 



dSidua , . 

— -— = 0, a = 1,2, 3; (38) 

OXj OXj 



In view of (|T8D we get from (|38|) the relations Six^ = 0, j = 1, 2, 3. Hence, 
-S*! = Si{t) and within the equivalence relation S (^) we may choose 5*1 to be 
equal to 0. Next, making use of the gauge invariance of SE (P we get from 
( |39D that within the equivalence relation S the equalities S2xi = 0, 2 = 1,2,3 
hold, whence 5*2 = S2{t). Again, in view of ( p!OD we may put 5*2 = 0, which 
means that the factor Q may be chosen as 1. With this result system (^) 
reduces to the system of three linear algebraic equations for the functions 

Al, ^2, ^3, 

i^at = -'^eUax.Ai, a = 1,2, 3. 

The determinant of this system is not equal to zero due to (|1^). Conse- 
quently, it has a unique solution. Making in this solution the hodograph 
transformation (^), we get the following expressions for Ai, A2, A^: 

^ 1 du{t,uj) 



2e dt 



Inserting into the above formula x = u{t, uj) from ( [2D| ) yields the explicit 
forms of the spacelike components of the vector potential of the electromag- 
netic field A{t,x), 

A{t, x) = — (^M (t) {x-w)+wj . (40) 

Here we use the designation 

M{t) = f{t)T-\t) + T{t)H{t)H-\t)T~\t), (41) 

where T{t), H{t) are variable 3x3 matrices defined by formulas ( PT| ) and 
(p2D, correspondingly, w = {wi(t),W2(t),W3(t))^ and the dot over a symbol 
means differentiation with respect to t. 
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Given the form ( ^OD of the spacehke components of the vector potential 
A{t,x), the condition for the magnetic field H not to vanish means that at 
least one of the expressions 

d + /?cos7, /3 cosa sin7 — 7 sin a, /3sina sin7 + 7 cosa (42) 

does not turn into zero. Hence, in view of the identity 

/ — (d + /3cos7) 

TT-^ = d + /3cos7 ^ 

\ f3 cos a sin 7 — 7 sin a /3 sin a sin 7 + 7 cos a 

— (/? cos a sin 7 — 7 sin a) 

— (/5 sin a sin7 + 7 cosa) 



whose validity is checked by straightforward computation, we conclude that 
T 7^ const. At last, solving (|l^) yields the explicit form of Aq, 

Ao{t, x) = -(j2 + To{t)] + lAa.^ - eAaAa. 

As the function Aq is real valued, we have to find the real part of the 
right-hand side of the above equality. Making use of (1^) yields 

3 hi 

'^e.Aaxa = X! I~- 

i=l "'i 

This, in its turn, gives the imaginary part of Tq = To(t), 

To = fo-'-j2^, Imfo = 0. (43) 



2^. hi' 

1=1 * 



Finally, we get 



1/3 f^i 1 f)( 1 \ 

Mt,x) = -(T. FUoo^)^^ + Ut)\ - eAaAa, (44) 
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where Fio(t^i), ^20(^^2), -^aolf^a) are arbitrary smooth functions, Ai, A2, A3 are 
given by (^OD, and the functions Ua = ujait, x), a = 1, 2, 3 are defined imphc- 
itly by formulas ([2D|)-(P^ and (^. In order to keep an exposition of the 



results compact, we do not present here Ao(t,x) in full detail, since the cor- 
responding expressions are too cumbersome. 
Thus we have proved the following assertion. 

Theorem 1 SE for the case of nonvanishing magnetic field admits sep- 
aration of variables if and only if the spacelike components Ai,A2, A3 of the 
vector potential of the electromagnetic field are linear in the spatial variables 
and given by and, furthermore, the timelike component Aq is given by 



Summing up, we conclude that conditions (^) and ( ^ID provide separabil- 
ity of SE for the case of a nonvanishing magnetic field. And what is more, the 
solutions with separated variables are of the form (|]) with Q = 1, where the 
functions uji{t, x),LJ2(t, x),uj3{t, x) are given implicitly by formulas (pOD-(p^) 
and (^). In fact, we have the eleven classes of vector potentials A{t, x) corre- 
sponding to the eleven classes of coordinate systems uja = uja(t, x), a = 1,2, 3. 
SE (d) for each class of the functions Ao{t, x), A{t, x) defined by (|0D and (||) 
under arbitrary Toit) , Fao{uJa) and fixed arbitrary functions a{t), [3{t), 7(t), 
Wa{i)i ha{t), a = 1, 2, 3 separates in exactly one coordinate system. The sep- 
aration equations read as (|11|)) where the coefficients Fai, a,i = 1, 2, 3 are the 
entries of the corresponding Stackel matrices (^0]), functions T^, a = 1,2,3 



are listed in (|3^ , and the functions Tq, Fao, a = 1, 2, 3 are arbitrary smooth 
functions defining the form of the timelike component of the vector potential 
A{t,x) (see (H). 

Note that proper specifying the functions Fao{uJa),ci = 1,2,3, may yield 
additional possibilities for variable separation in the corresponding SE. What 
we mean is that for some particular forms of the vector potential A{t, x) 



( ^41) there might exist several coordinate systems (POD-(P^) enabling us to 
separate the corresponding SE. However, the detailed study of this problem 
goes beyond the scope of the present paper. 

As an illustration of the previous subsection, we consider the problem of 
separation of variables in SE (|I]), where the vector potential of the electro- 
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magnetic field is of tlie form 



eA 



I -S2{t) \ 



Si{t) -S3(t) X, 



(45) 



V S2{t) Ssit) J 



eA, 







+ (si(t)a;i - S3{t)x3f + {S2{t)xi + S3{t)x2f) . 



Here q = const and 



2si(t) 

2S2(t) 

2s3(t) 



d{t) + Pit) cos7(t), 

cosa(t) sin7(t) — 7(t) sinQ;(t), 
I3{t) sinQ;(t) sin7(t) + 7(t) cosa(t), 



wliere Q;(t), 7(t) are arbitrary smootli functions. Evidently, choosing 
a(t)=const, /5(t)=const, and 7(t)=const yields the standard Coulomb po- 
tential. 

Making use of the results of Theorem 1, we conclude that SE with po- 
tential of the form ( ^5] ) separates in four coordinate systems, 

X = T(t)z, 

where T is the time-dependent 3x3 orthogonal matrix (|2T|) , and z is one of 
the following coordinate systems: 

1. spherical (formula 5 from (^), 

2. prolate spheroidal II (formula 6 from (^), where one should replace 
with ^3 = a(cothco'i tanhci;2 i l))? 

3. parabolic (formula 8 from (^), 

4. conical (formula 11 from (^). 

The separation equations ([TID for these cases take the form 
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II 



= -AiV?o, 

(A2 sech^u;2 - ^z)'^2, 



Integrating these equations yields a family of exact solutions of SE 
(|T|) with potential that are products of the exponential, confluent 
hypergeometric pO[) ^^^d Legendre pO] functions. 



ip'l = (XiO^ sinh~^ LUi — A2 sinh~^ cui - 
^2 = (Aia^ cosh"^ ct;2 + A2 cosh~^ uj2 

<^3 = ^3<^3- 



A3 + ga cosh tui sinh '^uji)(fi, 
- XsT QO. sinh UJ2 cosh~^ t^2)v^2, 



Integrating these equations yields a family of exact solutions of SE (|^) 
with potential (^51) in separated form that are products of the expo- 
nential and Coulomb spheroidal functions pi]. 



iipQ = -Aiv^o, 

= (Aie^'^1 - A2e2-i - A3 + 2ge2'^^)</.i, 
<^- = (Aie^'^^ + A2e2-2-A3)</^i, 
^3 = h^3- 

Integrating these equations yields a family of exact solutions of SE (H) 
with potential ( ^Sj) that are products of the exponential and confluent 
hypergeometric functions. 



= {\2k'^ Cn^{uJ2, k') - \-i)'^2, 

(fl = {X2k^ cn^(c^3, k) + A3)v?3- 
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Integrating these equations yields a family of exact solutions of SE 
([l|) with potential (|45|) that are products of the exponential, confluent 
hypergeometric, and Lame |]5D| functions. 



Note that these families of solutions are parametrized both by integration 
constants and by the three continuous spectral parameters Aj,i = 1,2,3. 
Under appropriate initial and boundary conditions the latter become discrete, 
i.e. Aj = Xin,n = 1,2,3,..., and we get a basis for expanding arbitrary 
smooth solutions of SE (|1]) with potential ( P5| ) in a properly chosen Hilbert 
space (for more details, see p^ ). 



2. The case of vanishing magnetic field. 



Provided condition {p7[j holds true, we can, without loss of generality. 



put A = at the expense of the gauge invariance of SE. Representing the 



complex-valued function Q in as Q = exp{Si+iS), inserting into (|T4D and 
splitting the obtained equations into real and imaginary parts analogously to 
what has been done for the case of a nonvanishing magnetic field we arrive 
at the conclusion that Si may be chosen to be equal to zero. Furthermore, 
the function S = S{t, x) satisfy the over-determined system of linear partial 
differential equations 

dS duJa dUa , . 

2t^7^ + ^ = 0, a = 1,2, 3. (46) 
oxj oxj at 

Solving it with respect to the derivatives dS/dxj, j = 1, 2, 3 (which is always 
possible due to (pTsP) and making in the relations obtained the hodograph 
transformation (^), yields 

2VS = M (t) {x-tv) + id, (47) 

where A^(t) is the 3x3 matrix (|4T|) and w = {wi(t) , W2it) , w^it))^ . The direct 
check shows that system ( ^Tj ) is compatible if and only if all the expressions 
given in (|42|) vanish identically. Hence, we conclude that the matrix T is 
constant. Utilizing the invariance of SE ([T|) with respect to the rotation 
group we can choose T = I, where / is the unit 3x3 matrix. Given this 
condition, (^) takes the form 

x = H{t)z{Lu)+w{t), (48) 
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and, furthermore, 

M{t) = H{t)H-\t). 

Next, integrating ( ^T]) and taking into account the equivalence relation 
(|TUp we get 

S-igg (I (49) 



Substituting (^91) into (|T5|) yields the form of AQ(t,x), 



the function T^it) being given by (^31) . 

Thus, we have proved the following assertion. 

Theorem 2 Given the restriction SE ^) admits a separation of vari- 
ables if and only if it is gauge equivalent to SE with ^ = and Aq being 
given by (|5Z 



Consequently, the conditions A = Q and (|50|) provide separability of SE 
for the case of vanishing magnetic field. Furthermore, the solutions with 
separated variables are of the form (|^) with Q = exp(iS'), where S = S(t,x) 
is given by (^U]), the functions u!i{t,x),LJ2(t,x),u!3(t,x) are given implic- 
itly by formulas (|48|), (0), and (H). Again, we have the eleven classes of 
vector potentials A{t,x) corresponding to the eleven classes of coordinate 
systems Ua = LJa(t,x), a = 1,2,3. SE (0) for each class of the functions 
AQ{t,x), A{t,x) = defined by (^) under arbitrary To{t), Fao{uJa) and fixed 
arbitrary functions Wa{t),ha{t), a = 1,2,3 separates in exactly one coordi- 
nate system. The separation equations read as (|ll]), where the coefficients 
Fai, a, i = 1,2,3 are the entries of the corresponding Stackel matrices (|30|), 
functions T^, a = 1, 2, 3 are listed in ( ^if ) and the functions Tq, F^q, a = 1, 2, 3 
are arbitrary smooth functions defining the form of the timelike component 



of the vector potential A{t,x) (see, (pO|)). 
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If we fix the temporal variable t to be equal to to ^ R- in the above ob- 
tained results, then it is not difficult to classify separable stationary Schro- 
dinger equations for a particle interacting with the electromagnetic field 

A{x) = {Ao{x),A{x)), 

{PaPa + eAo + E)^{x) = 0, (51) 

where Pa = id/dxa — eAa, a = 1, 2, 3, e = const, and £" is a spectral param- 
eter, and thus recover the classical result by Eisenhart [^. 

III. Separation of variables in the Hamilton— 
Jacobi equation 

It is well known that there exists a deep connection between the separation of 
variables in the Schrodinger and Hamilton-Jacobi equations (see, e.g., |I6| ). 
The Hamilton-Jacobi equation, 

ut + eAo + {u^^ + eAa) (u^, + eAa) = 0, (52) 

separates in any coordinate system providing separability of the Schrodinger 
equations (|I]) and, what is more, the inverse assertion is not true. We will 
make use of this connection for the sake of classifying separable Hamilton- 
Jacobi equations. 

First we fix the usual form of the separation Ansatz for the Hamilton- 
Jacobi equation, 

3 

u{t,x) = S{t,x) + ipo{t) + Lpi{uJi{t,x)), (53) 

i=l 

and, furthermore, fix the form of the ordinary differential equations for 

V50,V^l,V^2,V'3, 

ip'o = -Toit) - Ti{t)\i, if', = {-FaoM + Fa^MXif" . (54) 

Now, inserting the Ansatz (|53D into equation (^), eliminating the first 
derivatives of the functions yjo, V'l, V^2, V^s with the use of the above equa- 
tions and splitting by the variables ipQ, ipi, ip2, (ps, Xi, X2, ^s, we arrive at the 
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following system of nonlinear partial differential equations for the functions 



1^1^ = 0. .^j. ,, = 1.2.3; 

OXa OXa 

a =1.2.3; 

2||^ + e.')|..|. = 0. „=1.2.3; (55) 




j^-^ dXj dXj dt " dxa dxa dxa 

-To{t)+eAQ + ^AaAa = Q. 

The general solution uj = u{t, x) of the first six equations of the above 
system (which coincide with equations (|T2|) and (|13|)) can be reduced with 
the help of an appropriate equivalence transformation S to such a form that 
it satisfy the Laplace equation (||). 

It is not difficult to become convinced of the fact that making the change 
of variables, 

Q(t,x) = exp{iS{t,x)), (56) 

in (|T2])-(^) yields the system that coincides with (|^) with an exception 
of the last equation, where an additional term —i{AS + eAaxa) appears. As 
shown in the previous section, this term is a function of t only and is absorbed 
by Tq. Consequently, all the results on variable separation for SE apply to 
the case of the Hamilton-Jacobi equation ([5^) as well. 



IV. Concluding Remarks 



Theorems 1-2 give a complete solution of the problem of classification of 
SE's (P and (0) that are solvable within the framework of the method of 
separation of variables. By appropriate reductions of these results we can get 
the results on the separation of variables in SE's for a particle interacting with 
the electromagnetic field in one [H] and two ^ 



spatial dimensions. For 
example, in order to recover the results of |27| one has to choose dtp/dxs = 
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and consider the completely and partially split coordinate systems from the 
list i). 

It follows from Theorem 1 that the choice of magnetic fields H allowing 
for variable separation in the corresponding SE is very restricted. Namely, 
the magnetic field should be independent of spatial variables Xi,X2,x^ in 
order to provide the separability of SE (|l|) into three second-order ordinary 
differential equations. However, if we allow for separation equations to be of 
lower order, then additional possibilities for variable separation in SE arise. 
As an example we give the vector potential, 

A{t, x) = (^Ao ^xl + xi^ , 0, 0, As (^^xj + xl 

where Aq, A^ are arbitrary smooth functions. SE (|l|) with this vector-potenti- 



al separates in the cylindrical coordinate system t, ui = In yyxl + x^j , uj2 = 

arctan(xi/a;2), u;3 = x^ into two first-order and one second-order ordinary dif- 
ferential equations. The corresponding magnetic field H = rot A is evidently 
X dependent. 

As mentioned in the Introduction, a possibility of variable separation in 
SE is intimately connected to its symmetry properties. Namely, solutions 
with separated variables are common eigenfunctions of three mutually com- 
muting symmetry operators of SE. For all the cases of variable separation in 
SE ([l|) these operators can be constructed in explicit form, in analogy to what 
has been done in for the (l-|-2)-dimensional case. They are expressed 
in terms of the coefficients of the separation equations (|TT|). This fact en- 
ables application of the methods of the representation theory of Lie algebras 
for a further analysis of special functions arising as solutions of separation 



equations in the spirit of the famous Bateman's project [jl5 . 

The last remark is that the technique developed in the present paper 
can be directly applied in order to separate variables in the Pauli equation 
for a particle with spin | moving in the electromagnetic field and in the 
Fokker-Planck equation with a constant diagonal diffusion matrix. 

A study of the above-mentioned problems is in progress now and will be 
reported in our future publications. 
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